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ABSTRACT 

We present a spectroscopic analysis of nearly 8000 late-type dwarfs in the Sloan Digital Sky Survey. Using the 
Ha emission line as an activity indicator, we investigate the fraction of active stars as a function of spectral type 
and find a peak near type M8, confirming previous results. In contrast to past findings, we find that not all M7-M8 
stars are active. We show that this may be a selection effect of the distance distributions of previous samples, as 
the active stars appear to be concentrated near the Galactic Plane. We also examine the activity strength (ratio of 
the luminosity emitted in Ha to the bolometric luminosity) for each star, and find that the mean activity strength 
is constant over the range M0-M5 and declines at later types. The decline begins at a slightly earlier spectral 
type than previously found. We explore the effect that activity has on the broadband photometric colors and 
find no significant differences between active and inactive stars. We also carry out a search for subdwarfs using 
spectroscopic metallicity indicators, and find 60 subdwarf candidates. Several of these candidates are near the 
extreme subdwarf boundary. The spectroscopic subdwarf candidates are redder by ~ 0.2 magnitudes in g — r 
compared to disk dwarfs at the same r—i color. 

Subject headings: solar neighborhood — stars: low-mass, brown dwarfs — stars: activity — stars: late-type — 
subdwarfs — Galaxy: structure 



1. INTRODUCTION 

Low mass stars with late spectral types (M,L) are the major- 
ity constituent of the Galaxy by number. Their main-sequence 
lifetimes are much greater than the current age of the Universe 
and they therefore serve as useful probes of Galactic star for- 
mation history in the local solar neighborhood (Gizis, Reid & 
Hawley 2002). They also encompass many important regions of 
stellar parameter space, including the onset of complete convec- 
tion in the stellar interior, the onset of significant electron de- 
generacy in the core, and the formation of dust and subsequent 
depletion of metals onto dust grains in the stellar atmosphere. 
Of particular interest is the fact that many late-type stars have 
strong surface magnetic fields (Johns-Krull & Valenti 1996) 
that heat the outer atmosphere above the photosphere, and lead 
to observable emission from the chromosphere (e.g. Ca II and 
H Balmer series lines), the transition region (e.g resonance lines 
of abundant ions such as C IV), and the corona (e.g. thermal 
soft X-rays). The physics that controls the production of mag- 
netic fields in low mass stars is not well understood, as the lack 
of a radiative-convective boundary layer precludes storing large 
scale fields as in a solar-type dynamo. However, recent work on 
turbulent dynamo mechanisms (Bercik et al. 2004) may soon 
provide new insight into the formation and properties of mag- 
netic fields in low mass stars. The present study seeks to place 
strong empirical constraints on the magnetic activity in these 
stars, by measuring optical Ha emission in a sample of nearly 
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8000 low mass dwarfs. Our results will aid in the interpre- 
tation of the models and provide a connection to the physical 
processes occurring in the atmospheres of the stars. 

The Ha emission line is produced by collisional excitation 
in the relatively dense chromospheres of these low mass, high 
gravity stars. Located in the red region of the optical spectrum, 
it is the strongest and best-studied indicator of magnetic activ- 
ity in late type stars (in contrast to solar type stars which are 
usually studied in the Ca II H and K resonance lines). Previous 
results indicate that the fraction of M dwarfs with Ha emission 
increases from early to mid M spectral types (Joy & Abt 1974; 
Hawley, Gizis & Reid 1996), reaches a peak near type M7 
where essentially all stars are active, and declines toward later 
types (Gizis et al. 2000a). The increasing fraction in the early 
to mid M types may be an age effect such that the activity lasts 
longer in the mid M types (Hawley, Reid & Tourtellot 2000a; 
Gizis et al. 2002), while the decline toward late M types may 
reflect difficulty in producing and/or maintaining surface mag- 
netic fields due to the physics of turbulent dynamos and/or the 
increasingly neutral atmosphere (Hawley, Reid & Gizis 2000b; 
Fleming, Giampapa & Schmitt 2000; Mohanty et al. 2002). It 
does not appear that the atmosphere changes character such that 
Ha is no longer produced when magnetic activity is present in 
the later type stars. Indeed the ratio of Ha to transition region 
(C IV) and soft X-ray emission appears to be similar throughout 
the M dwarf spectral sequence (Hawley & Johns-Krull 2003). 

Previous work also shows that the activity strength, measured 
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by the ratio of the luminosity in Ha to the bolometric luminos- 
ity (L Ha /L/, ;), is nearly constant (with large scatter) through the 
M0-M6 range (Hawley et al. 1996) and declines at later types 
(Burgasser et al. 2002, Cruz & Reid 2002). However, the sam- 
ple of stars on which these conclusions are based numbers less 
than 100 for types M7 and later, due to the difficulty of obtain- 
ing spectra for these faint objects. Our sample in this crucial 
spectral type range is an order of magnitude larger, compris- 
ing more than 1000 stars with types between M7 and L0, and 
nearly 8000 stars between M0 and L0. This uniformly acquired 
and reduced sample from the Sloan Digital Sky Survey (SDSS), 
with well-characterized uncertainties, now allows us to obtain 
statistically significant magnetic activity results for the entire 
M dwarf spectral type range. In addition, we compare the pho- 
tometric colors of the active and inactive stars to investigate 
whether the active stars are typically bluer at the same spectral 
type, an effect attributed to an increase in plage areas on active 
stars by Amado & Byrne (1997). Alternatively, such a corre- 
lation might indicate that micro-flaring activity plays a signif- 
icant role in heating the outer atmospheres of active M dwarfs 
(cf. Gudel et al. 2003). 

Our SDSS spectral sample is also well-suited for investiga- 
tion of the incidence of low mass subdwarfs in the local neigh- 
borhood. Gizis (1997) showed that M subdwarfs could be iden- 
tified by comparing the metallicity and gravity sensitive molec- 
ular bands of CaH and TiO. Using his technique, we identify 
dozens of new M subdwarf and extreme subdwarf candidates, 
and examine their colors in the SDSS ugriz filters. 

2. DATA 

The SDSS (York et al. 2000; Gunn et al. 1998; Fukugita et 
al. 1996; Hogg et al. 2001; Smith et al. 2002; Stoughton et 
al. 2002; Pier et al. 2003, Abazajian et al. 2003) is an excel- 
lent tool for spectroscopic studies of stars. Each spectroscopic 
plate produces 640 spectra with resolution A/ A A = 1800. Be- 
cause all spectra are obtained and reduced through uniform 
pipelines with well understood uncertainties (York et al. 2000), 
it is possible to undertake statistically robust studies of very 
large samples. Our sample was selected from the SDSS spec- 
troscopic and photometric databases as follows. In July 2002, 
the spectroscopic database was queried for all objects which 
satisfied the SDSS color cuts for late-type stars given in Haw- 
ley et al. (2002; hereafter H02), viz. 0.8 < (r-i) < 3.0 and 
0.3 < (i — z) < 3.5. The photometry for the color cuts came 
from the preliminary SDSS photometric database 9 available at 
that time. This color cut yielded ~ 17000 late-type star candi- 
dates with SDSS spectra. In April 2003 we added an additional 
~ 5000 candidates from spectroscopic observations that were 
specifically targeted for late type stars in the Southern Survey 
region. The final ~ 22000 candidates come from 427 SDSS 
spectroscopic plates; 380 of the plates are in the SDSS Data Re- 
lease 2 (DR2; Abazajian et al. 2004) sky area while 47 plates 
have yet to be publicly released. Due to recent improvements 
in the photometric reductions, we have used photometry only 
from the 5_4_25 version of Photo (SDSS photometric reduc- 
tion software) in our final analysis. This is the same version 
used in the DR2 release. 

3. ANALYSIS 
3.1. Spectral Types 

9 the "Chunk" database 



In order to assign spectral types to the stars in the sample, 
we pass each spectrum through two independent spectral typ- 
ing pipelines. The first method determines the best fit, in a 
least squares sense, to both the SDSS late-type template spec- 
tra defined in H02 and the template spectra from the Keck 
and Apache Point Observatories used by H02 in their analy- 
sis; it is a derivative of the procedure originally described in 
Kirkpatrick (1991) and subsequently modified by Henry et al. 
(2002). The second method measures several molecular band 
indices in each spectrum and uses a weighted mean of the in- 
dividual types from each index to assign a best-guess spectral 
type (cf. Reid, Hawley, & Gizis 1995a; Kirkpatrick et al. 1999). 
Both methods are described in detail in H02. Using the spectral 
typing pipeline outputs as a guide, each spectrum is individu- 
ally inspected by eye and the assigned spectral type is altered 
as needed. As reported in H02, the uncertainty in this proce- 
dure is <~ 1 spectral type. Spectra with a signal-to-noise ratio 
(SNR) < 5 are not typed, as experience shows that they do not 
give useful results in our analysis. Of the 22000 original can- 
didates, ~ 9200 were assigned M or L spectral types and had 
spectra with the necessary signal-to-noise ratio for further anal- 
ysis. An important byproduct of the spectral typing procedure 
was the identification of numerous white dwarf - M dwarf bi- 
nary systems in our sample. We find the binaries based on the 
appearance of blue continuum flux and features such as broad 
hydrogen Balmer absorption lines (Raymond et al. 2003). We 
have removed all such binaries from the final sample so that 
their colors will not affect our investigation of the colors of ac- 
tive vs. inactive M dwarfs, and because the nature of the M 
dwarf activity may be influenced by the presence of a close bi- 
nary companion. §3.3 outlines a photometric test used to further 
constrain our sample to stars without white dwarf companions. 
Our sample is thus, as far as possible, composed of only single 
M dwarfs (or M dwarfs in binary systems with fainter M or L 
dwarfs that we are unable to discern from the spectra). 

3.2. Activity 

3.2. 1 . Ha Measurements 

The SDSS spectra are spectrophotometrically calibrated, so 
we measure the Ha flux directly from each spectrum. The chal- 
lenge is to automatically, but accurately, measure the flux in Ha 
over a range of spectral types with changing continuum slope 
and molecular features. Our procedure uses trapezoidal inte- 
gration to sum the flux under the emission line, taking care to 
avoid and/or account for regions that are affected by local TiO 
features including the continuum bump just to the blue of Ha 
(Pettersen & Coleman 1981) and the absorption trough to the 
red (Tinney & Reid 1998; White & Basri 2003). We define a 
14A interval (7 A on either side of the line center) for the Ha 
flux calculation. We investigated several options to find the op- 
timum continuum flux measure: a) the average over two 50A 
sections (6500-6550A and 6575-6625A) on either side of the 
line; b) the average over two 10A sections (6545-6555A and 
6575-6585A) on either side of the line; and c) the blue 10A 
region only, as used by Reid et al. (1995a) for the PMSU sur- 
vey (6545-6555A). We found no significant differences in the 
mean continuum values for these three measures, and therefore 
chose to use measure (a) as it is the most robust, and because 
we also use the 50A regions as comparison for the equivalent 
width (EW) measurement (see below). The formal uncertainty 
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in the flux at each wavelength is available from the SDSS spec- 
tral calibration and is propagated at every stage of the analysis. 

In order to be defined as active, each spectrum must pass the 
following criteria: 1) the Ha EW must be greater than 1.0A; 2) 
the uncertainty in the EW must be smaller than the EW value; 
3) the height of the Ha emission line must be 3 times the noise 
at line center; and 4) the EW of Ha must be 3 times larger than 
the average EW of the 50A comparison regions. We exam- 
ined subsets of the spectra that pass less conservative criteria, 
including those that passed only 2 or 3 of the 4 tests. These ac- 
count for less than 5% of the total sample, and upon inspection 
are clearly marginal. There are 430 stars in the less conserva- 
tive subsets that appear to be active but have Ha EW less than 
1 .OA. We do not include these in our analysis but mention them 
as potential active stars similar to the stars in Groups C and E 
of Gizis et al. (2002). For the remainder of this paper all results 
reported will be from spectra that either passed all of the above 
criteria (active - 1910 stars) or failed all of our tests, including 
the tests with less stringent criteria (inactive - 5930 stars), for a 
total of 7840 stars in the final activity sample. Subdwarfs have 
also been removed from this sample; they are discussed sepa- 
rately in §4.4 below. We reiterate that in order to be classified 
as inactive, the Ha EW must be less than 1 .OA. Thus a star with 
Ha EW > 1 .OA, but with noisy data, is not classified as either 
active or inactive, but is simply removed from the sample. The 
H02 subset of 659 late type stars from the SDSS Early Data 
Release was analyzed with this method and then inspected by 
eye. We found that less than 4% of the stars were incorrectly 
labeled as active or inactive, and cite this as a measure of our 
ability to determine the activity status automatically. 

3.2.2. Quantifying Activity 

A star's magnetic activity can be qualitatively assessed as ei- 
ther active or inactive using the EW of the Ha emission line. 
However, because the EW measure depends on the flux in 
the local stellar continuum, it cannot be used to compare the 
strength of the activity in stars of different spectral type (Reid, 
Hawley, & Mateo 1995b). Instead, we measure activity strength 
by computing the ratio of the luminosity in Ha (L Ha ) to the 
bolometric luminosity (L/, ;)- The ratio Ln a fLk [ is a measure 
of the magnetically driven radiative output compared to the to- 
tal radiative output of the star, and does not depend on the shape 
of the stellar continuum. 

In the past, L Hq has been calculated by converting the mea- 
sured Ha flux to a luminosity using the distance to the star. 
However, distance determinations for the stars in our sample 
rest on either spectroscopic or photometric parallax methods. 
To avoid introducing additional uncertainty into our L Hq calcu- 
lation due to uncertainties in the distances, we use a distance- 
independent method of calculating Lh q /L/, /. The method is 
described in detail in Walkowicz et al. (2004) and uses well- 
observed nearby M and L dwarfs with known distances to em- 
pirically determine the ratio of the luminosity in the continuum 
at Ha, L(cont)H Q , to the bolometric luminosity. The resulting 
ratios \ = L(cont)H Q /Lfe ; are tabulated as a function of spec- 
tral type and color, and when multiplied by the EW in Ha for a 
given star, give the desired ratio L Ha /L/, ;. 

Since the \ ratios are tabulated at 0.5 spectral type intervals, 
and we use integer spectral types in our analysis, a weighted 
mean of the \ values surrounding a given integer spectral type 
was computed. We give the two surrounding values 1/4 the 
weight of the integer \ value. This weight is also used in prop- 



agating the uncertainties. We multiply each measured Ha EW 
by the appropriate mean x factor and thereby calculate Lu a fLb i 
for each active star in our sample. Note that because there are 
no M0 stars in the nearby star sample, the M0 Lu a /^boi values 
in the present work have been calculated using only the M0. 5 \ 
value. 

3.3. Colors 

The colors used in our analysis are formed from the ugriz psf 
magnitudes in the SDSS photometric database. As described 
above, we use photometry only from the latest (Fall 2003, Photo 
v. 5_4_25) pipeline reductions. To assure reasonable photomet- 
ric accuracy, we include only stars with u magnitude uncertain- 
ties less than 0.3 and griz magnitude uncertainties less than 0.2 
in our analysis of color effects. This decreases the number of 
stars available for analysis to ^2200. However, if the u band 
data are not required, for example when examining (g — r) col- 
ors, the sample numbers ^5800. The maximum uncertainties in 
the colors that we use are therefore 0.36 magnitudes in (u-g) 
and 0.28 magnitudes in {g — r), (r — i), and (i — z). The great 
majority of stars in this color sample have uncertainties much 
smaller than these maximum values, with typical values being 
a(r) < 0.02 for r < 17 (Ivezic et al. 2003; 11% of the sample); 
a(r) < 0.03 for r < 19 (46% of the sample); and a(r) < 0.04 
for r < 21 (83% of the sample). 

As discussed above, a significant source of scatter in the blue 
colors of the late type stars in our sample was introduced by 
the presence of numerous white dwarf - M dwarf binary sys- 
tems. The spectroscopically-identified pairs were originally re- 
moved during the spectral typing process. In addition, we made 
a color cut on our sample, and included only stars whose col- 
ors lie outside the locus of known white dwarf - M dwarf pairs 
in the (g — r) vs. (u-g) color-color diagram (Silvestri, Hawley 
& Szkody 2003; Smolcic et al. 2004). In practice this meant 
that only stars with (u-g) > 1.8 are included in the sample. 
All suspected pairs based on either photometry or spectroscopy 
were therefore removed, and are not included in the final spec- 
troscopic sample of 7840 stars. 

3.4. Molecular Band Indices 

To investigate the presence of metal poor stars in our sam- 
ple, we use the varying strengths of the TiO and CaH molecular 
bands in a low metallicity environment. We measure the Ti05, 
CaHl and CaH2 molecular band indices following the proce- 
dure outlined by Reid et al. (1995a). Gizis (1997) uses the same 
indices to define divisions between dwarfs, subdwarfs and ex- 
treme subdwarfs. In our study, the following criteria are used to 
restrict the sample of subdwarfs to high probability candidates: 
a) the band indices must be larger than their uncertainties (typ- 
ical uncertainties are 0.04 in Ti05; 0.03 in CaHl; and 0.02 in 
CaH2); b) the band heights must be larger than both the noise 
in the continuum near the bandhead and the noise in the bottom 
of the absorption feature; and c) the computed uncertainties in 
the band indices must not overlap the subdwarf or extreme sub- 
dwarf boundaries. 

4. RESULTS 

4. 1 . Activity Fraction vs. Spectral Type 

We find that in our sample of 7840 stars that are either ac- 
tive or inactive, 1910 are classified as active (24.4%). Figure 1 
shows the fraction of stars which are active as a function of 
spectral type. The active fraction peaks at spectral type M8, 
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where 73% of stars are active. The numbers above each spectral 
type bin represent the total number of stars (both active and in- 
active) in the bin. The peak at M8 is consistent with Gizis et al. 
(2000a), who found a peak in the relation at types M7 and M8 
using similar activity criteria. However, all stars in their sample 
were active at these types, while our results show ^64-73% 
of our M7-8 stars are active. We have confirmed by eye that 
indeed many of the M7 and M8 stars in our sample do not have 
measurable activity. Figure 2 shows examples of active and in- 
active M7 and M8 stars. The continuum features near Ha are 
noticeable in the inactive stars (blue bump, red trough), but it is 
clear that no measurable emission at Ha is present. 

Following a suggestion by Neill Reid (private communica- 
tion), we considered whether the disparity in active fraction at 
type M7 could be a distance effect. The Gizis et al. (2000a) 
sample was drawn from 2MASS and represents mostly nearby 
stars, while the SDSS sample described here extends to much 
greater distances due to the fainter limiting magnitude. We 
therefore restricted ourselves to just the M7 stars in the sample, 
and computed the distance to each star via the (i—z) photomet- 
ric parallax relation from H02. Only stars with good photom- 
etry (a(i-z) < 0.28) were included (see §3.3), which reduced 
the sample from 777 to 269. Note that we did not combine 
stars of different spectral type (color) in the sample, so the rela- 
tive distances between the M7 stars should be quite accurate, as 
uncertainties in the absolute calibration of the photometric par- 
allax relation at different colors do not affect the calculation. 
We then computed the vertical distance, perpendicular to the 
Galactic Plane, for each star. The results are shown in Figure 3 
for 60 pc distance bins. Though the number of stars in each bin 
is relatively small, leading to significant Poisson errors, there 
does appear to be a trend of decreasing active fraction with ver- 
tical distance from the Plane in both the North (positive) and 
South (negative) directions. If activity is related to age, as dis- 
cussed by Hawley et al. (2000a) and Gizis et al. (2002), then 
this effect would be a natural consequence of the increasing ve- 
locity dispersion in older disk stars: the older stars would be 
more likely to be found further from the Plane where they were 
born. Kirkpatrick et al. (1994) first discussed the possibility 
that younger stars were found when looking South (toward the 
Plane, from the Sun's position 30 parsecs North of the Plane) 
than North, out of the Plane. Small sample sizes and uncer- 
tain distances precluded a definitive test using data available at 
that time, but we may now be seeing evidence for this effect in 
our much larger sample. The uncertainties are still large how- 
ever, and more data are needed to verify this result. We plan to 
carry out a rigorous comparison between the 2MASS and SDSS 
samples to further investigate the correlation of activity fraction 
with distance from the Plane. 

4.2. Activity Strength vs. Spectral Type 

Figure 4 illustrates the activity strength, log(L HQ /Lfoo/X as a 
function of spectral type. Each histogram represents the distri- 
bution of log(L HQ /L/, ;) at a given spectral type. The number 
to the right of the histogram is the peak value of the distribu- 
tion. We see that log(L Ha /Li ;) peaks near -3.7 for stars of 
type M0-M5, in agreement with Hawley et al. (1996), and sys- 
tematically peaks at lower activity strength among stars with 
later type. The 1A Ha EW measurement limit in our sample 
imposes a lower limit on the measurable activity strength. At 
each spectral type the last bin (smallest log(L Ha /L/, ;)) with 
stars in it corresponds to this measurement limit. Table 1 gives 



the log(LH Q /Lfc /)) corresponding to an Ha EW of 1 A for each 
spectral type. Thus, the broader distributions seen at later types 
may be partly due to a selection effect. However, the distribu- 
tions among the later types are quite symmetric and definitely 
declining at log(LH Q /L/, ;) values that are easily measurable. 
Therefore it appears that the activity strength covers a larger 
range, and is less sharply peaked, among the later type stars. 

Figure 5 shows the mean of the L Ha /Lboi distribution as a 
function of spectral type. The errors in this plot are the errors 
in the Lua/^boi calculation added in quadrature with the la- 
distribution of each bin. The dotted line indicates the L Hq /L/, / 
for an Ha EW of 1A (the required minimum EW from §3.2). 
Table 1 gives the mean values for each spectral type. The ap- 
proximately constant value from M0-M5 is clearly shown, as is 
the decline starting at types M5-M6. This decline occurs at a 
slightly earlier type than in previous studies, where the decline 
appeared at types M7-M8 (Gizis et al. 2000a; Burgasser et al. 
2002). 

4.3. Colors and Activity 

We used several SDSS color-color diagrams to investigate 
the photometric properties of active stars. Figure 6 shows that 
the colors of the active stars are very similar to the inactive stars. 
It has been claimed previously that active stars are marginally 
bluer than inactive ones in the Johnson U-B color, with <5 (U- 
B) = 0.04 ± 0.09 for late-type main-sequence dwarfs (Amado 
& Byrne 1997). They speculate that the effect could be due 
to the appearance of bright plage regions on the active stars. 
Our results for the SDSS (u-g) color show a similar, but not 
statistically significant effect: S(u-g) = 0.05 ± 0.25. In the 
(g - r) color the effect has the opposite sign: 5(g - r) = -0. 1 1 ± 
0.16. The scatter in the color distribution is much larger than 
the mean offsets, and we therefore cannot confirm the hypoth- 
esis that plage regions are more prevalent on the active stars. A 
two-dimensional Kolmogorov-Smirnov test confirms that there 
is no significant difference between the active and inactive stars 
in any of the color-color diagrams. 

4.4. Subdwarfs 

Figure 7 (left panel) plots the CaHl index versus the Ti05 
index. The solid line denotes the dwarf/subdwarf boundary de- 
fined by Gizis (1997). The CaHl-Ti05 relation is valid only for 
Ti05 > 0.49 and therefore we have not extended the boundary 
beyond that value. We find 60 stars (crosses) that fall below the 
subdwarf boundary (and have uncertainties that do not cross the 
boundary as described above). Figure 7 (right panel) shows a 
similar plot for the CaH2 index. Again the solid line indicates 
the dwarf/subdwarf boundary and the dashed line is the bound- 
ary between the subdwarf and extreme subdwarf populations. 
Several of the stars are close to the extreme subdwarf bound- 
ary and warrant further observation. The subdwarf and extreme 
subdwarf candidates from this study are listed in Table 2. 

The colors of the subdwarf candidates are compared to the 
normal disk population in Figure 8. The subdwarfs are clearly 
delineated in the (g-r) color, with the median (g — r)=1.61, 
nearly 0.2 magnitudes redder than the normal disk dwarfs. This 
effect may be due to excess absorption in numerous hydride 
bands in the g bandpass, and has been seen previously in the 
Johnson (B-V) color (Hartwick 1977; Bergbusch & Vanden- 
Berg 1992; Gizis 1997; Dahn et al. 1995). Figure 9 shows the 
strong correlation between (g-r) color and CaHl band strength 
for stars with r < 19.5. The subdwarf candidates are clearly 
separated from the normal disk dwarfs in this diagram. This 
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allows subdwarf candidates to be chosen photometrically, for 
example in the SDSS spectral targeting pipeline, with reason- 
able success. Note that not all of the stars found here come 
from subdwarf targeting; some were serendipitous observations 
of objects targeted for other purposes. 

5. SUMMARY 

Using a sample of more than 22000 candidates from the 
SDSS spectroscopic database, 7840 late type stars were iden- 
tified and used to investigate the magnetic activity in M and 
early L dwarfs. This sample is much larger than all previous 
samples used for this purpose. Our results show that: 

1. The fraction of active stars rises monotonically from 
spectral type MO to M8 and then declines monotonically 
to the latest types we measured (L3-L4). 

2. Only ~ 70% of the M7 stars in our sample are active, 
which differs from previous studies showing 100% activ- 
ity at spectral type M7. We suggest that the answer to 
this discrepancy may lie in the distance distributions of 
the samples used. Our results indicate the possibility that 
stars near the Galactic Plane are more likely to be active, 
as expected if activity depends on age. 

3. The activity strength, as measured by Lh q /L/, ;, is ap- 
proximately constant from types M0-M5, and declines at 
later types. The width of the distribution is relatively nar- 
row at early types but broadens considerably at type M6 
and later. 

4. The average (u-g) color is slightly bluer for active stars 
compared to inactive ones, while the average (g—r) color 
is slightly redder. However, these differences are not sta- 
tistically significant. 

5. 60 new subdwarf candidates are identified spectroscop- 
ically by comparing the strengths of the TiO and CaH 
molecular bands. Several of the candidates lie near the 
extreme subdwarf boundary and may be interesting tar- 
gets for future observations. 

6. The subdwarf candidates have similar colors to the disk 
dwarfs except in the (g - r) color, where the median 
subdwarf (g - f) is almost 0.2 magnitudes redder than 
disk dwarfs at the same (r—i) color. This effect allows 
subdwarfs to be chosen photometrically from the SDSS 
database. 
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Fig. 1. — The fraction of active stars is shown as a function of spectral type. The numbers above each bin represent the total 
number of stars used to compute the fraction. The peak at spectral type M8 is consistent with previous results, but the presence of 
inactive M7 and M8 stars was not seen in earlier observations. Spectral type bins with no observed activity (L3, L4) are included for 
completeness. 
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Fig. 2. — Active and inactive M7 (left panels) and M8 (right panels) stars. Small boxes in each panel contain a 100A section of the 
spectrum near Ha. The dashed lines are drawn at the wavelength of Ha. Note that the line falls to the right of the closest spectral 
feature in both inactive spectra. 
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Fig. 3. — The fraction of M7 stars that are active is shown as a function of vertical distance above and below the Galactic Plane. 
Active stars appear to be concentrated toward the Plane, which may be an age effect (see text). 
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Fig. 4. — The distribution of the activity strength, log(LHa/L/*>/), is shown as a function of spectral type. The numbers to the right 
at each spectral type give the peak value of the histogram. The mean activity strength clearly declines at types M6 and later, and the 
distribution appears to broaden. 
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Fig. 5. — The mean activity strength, log(mean L Hct /L fco /), is plotted as a function of spectral type. Uncertainties reflect the propa- 
gated errors in the calculation of Lh q /1-boi added in quadrature with the 1 a distribution of Lh q /L/ m ; in each bin. The decline in mean 
activity strength begins at spectral types M5-M6. The dotted line indicates the Lua/^boi value that we are sensitive to, due to the 1 A 
EW activity criterion. 
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Fig. 6. — Several color-color diagrams in the SDSS ugriz filters are used to compare the colors of active stars (diamonds) to inactive 
stars (dots). A 2D KS test shows that there are no significant differences between active and inactive stars in any of the colors. 
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Fig. 7. — The left panel compares the CaHl and Ti05 molecular band strengths, while the right panel compares CaH2 and Ti05. The 
solid line in each panel is the dwarf/subdwarf boundary from Gizis (1997). Crosses are subdwarf candidates that meet the criteria 
defined in S3. 4. 
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Fig. 8. — These SDSS color-color diagrams compare the colors of subdwarf candidates (diamonds) with normal disk dwarfs (dots). 
The (g-r) colors of the subdwarf candidates are systematically redder at a given (r-i) and (i-z) color. 
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Fig. 9. — The redder (g — r) colors of the subdwarf candidates are clearly evident in the (g — r) vs. CaHl diagram. The SDSS (g — r) 
color provides a good photometric proxy for measuring strong hydride bands and allows subdwarfs to be identified photometrically. 
To reduce the photometric scatter, only stars with r < 19.5 are shown. 



Table 1 
Lna/Lboi vs. Spectral Type 



Spectral Class 


log(meanL Hc( /L 6o ;) 


log(L Hct /U ,) for lA EW 


MO 


-3.67 


-3.94 


Ml 


-3.62 


-3.93 


M2 


-3.63 


-4.02 


M3 


-3.60 


-4.13 


M4 


-3.57 


-4.20 


M5 


-3.67 


-4.56 


M6 


-4.00 


-4.76 


M7 


-4.31 


-5.28 


M8 


-4.10 


-5.22 


M9 


-4.22 


-5.42 


LO 


-4.37 


-5.33 
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Table 2 
Subdwarf Data 



NT 

Name 




,„ „ab 

g~ r 


;ab 

r— i 


l — Z 




CaHl 


CaHz 


Ha hSW ( A) 


Sp. lype 


SDSS J002228.00-091444.8 


18.92 


-0.89 


3.44 


0.48 


0.65 


0.71 


0.52 


1.16 


1 


SDSS J00354 1.84+0032 10.1 


20.15 


1.64 


1.14 


0.59 


0.68 


0.69 


0.54 


-0.16 


3 


SDSS J003701. 37-003248.3 


20.15 


1.60 


1.23 


0.67 


0.57 


0.64 


0.46 


-0.23 


3 


SDSS J003755.20-002134.2 


18.70 




1.01 


0.89 


0.65 


0.75 


0.57 


-0.35 


1 


SDSS JO 108 11.89+003042.4 


17.34 


1.61 


1.17 


0.65 


0.51 


0.71 


0.45 


-0.43 


3 


SDSS J012853. 17+003356.6 


19.34 


1.56 


1.12 


0.58 


0.58 


0.69 


0.50 


-0.21 


3 


SDSS J014631. 67+001658. 3 


19.79 


1.26 


1.39 


0.49 


0.58 


0.72 


0.48 


-0.18 


3 


SDSS J020753.57-001958.8 


20.16 


1.54 


1.22 


0.64 


0.50 


0.67 


0.46 


0.02 


3 


SDSS J021956.49+005153.6 


19.27 


1.60 


0.97 


0.45 


0.72 


0.74 


0.54 


2.97 


2 


SDSS J022221. 54+005430.3 


19.01 


1.64 


1.13 


0.59 


0.56 


0.70 


0.46 


0.04 


3 


SDSS J024501. 77+0033 15.8 


19.38 


1.54 


1.03 


0.58 


0.64 


0.73 


0.51 




1 


SDSS J031314.28-000619.8 


20.33 


1.48 


1.22 


0.69 


0.51 


0.66 


0.51 


-3.27 


2 


SDSS J032749.80-004445.4 


19.88 


1.59 


0.99 


0.55 


0.64 


0.70 


0.62 


-0.35 


2 


SDSS J033408.64-072349.2 


20.21 


1.89 


0.91 


0.55 


0.79 


0.65 


0.45 


-0.85 


2 


SDSS J081329.95+443945.6 


19.39 


1.68 


1.07 


0.60 


0.63 


0.70 


0.50 


-0.45 


2 


SDSS J082230.00+471645.8 


19.39 


1.67 


1.15 


0.56 


0.54 


0.71 


0.46 


0.04 


2 


SDSS J083002.73+483251.6 


19.92 


1.90 


1.17 


0.62 


0.53 


0.69 


0.45 


1.89 


3 


SDSS J083217.77+522408.2 


19.61 


1.49 


1.10 


0.55 


0.60 


0.71 


0.52 


-0.38 


2 


SDSS J084105.39+032109.6 


20.05 


1.45 


1.15 


0.58 


0.62 


0.69 


0.47 


0.24 


3 


SDSS J085843.89+51 1210.1 


20.00 


1.53 


1.05 


0.52 


0.66 


0.70 


0.60 


-0.98 


2 


SDSS J090238.75+471813.6 


19.97 


1.73 


1.02 


0.46 


0.66 


0.71 


0.64 


-0.58 


2 


SDSS J090434.02+513153.9 


19.40 


1.78 


0.93 


0.58 


0.78 


0.61 


0.45 


-0.64 


1 


SDSS J091451.98+453152.8 


19.02 


1.64 


1.23 


0.62 


0.51 


0.70 


0.50 


-0.58 


3 


SDSS J092429.76+523410.7 


18.70 


1.67 


0.99 


0.55 


0.61 


0.63 


0.45 


-0.47 


2 


SDSS J092534. 16+524442.4 


19.79 


1.71 


0.95 


0.65 


0.71 


0.73 


0.58 


-0.61 


3 


SDSS J092708. 10+561648.1 


19.37 


1.62 


0.93 


0.47 


0.71 


0.70 


0.62 


-0.41 


1 


SDSS J092745.78+582122.7 


20.47 


1.78 


1.21 


0.57 


0.59 


0.58 


0.55 


-0.86 


3 


SDSS J093024.66+554447.7 


19.14 


1.56 


1.13 


0.27 


0.57 


0.70 


0.48 


-0.21 


3 


SDSS J093141.85+453914.5 


19.45 


1.58 


0.96 


0.54 


0.61 


0.69 


0.57 


-0.41 


2 


SDSS J094306.37+465701.4 


19.74 


1.59 


1.06 


0.64 


0.54 


0.61 


0.50 


-1.34 


3 


SDSS J095147.77+003612.0 


18.27 


1.58 


1.04 


0.53 


0.62 


0.74 


0.53 


-0.35 


2 


SDSS J100109.54+015450.2 


19.12 


1.71 


0.98 


0.56 


0.66 


0.72 


0.45 


-0.62 


2 


SDSS J101031. 13+651327.6 


19.39 


1.68 


1.05 


0.57 


0.68 


0.69 


0.44 


2.80 


1 


SDSS J104320.47+010439.4 


19.16 


1.66 


1.11 


0.60 


0.70 


0.66 


0.44 




2 


SDSS J105122.43+603844.8 


17.15 


1.65 


1.00 


0.56 


0.67 


0.62 


0.44 


-0.17 


2 


SDSS Jl 12751.35-001246.8 


20.03 




1.07 


0.65 


0.66 


0.73 


0.56 




1 


SDSS Jl 13501.76+033720.3 











0.53 


0.64 


0.53 





3 


SDSS Jl 15900.70+665214.3 


19.37 


1.61 


1.09 


0.65 


0.58 


0.68 


0.47 


-0.53 


3 


SDSS J120724.85+004346.4 


18.10 


1.49 


1.19 


0.62 


0.52 


0.72 


0.49 


-0.44 


3 


SDSS J121510.41+003342.7 


18.00 


1.41 


0.70 


0.45 


0.87 


0.84 


0.77 


-0.19 





SDSS J125919.29-025402.3 











0.60 


0.74 


0.52 


-0.39 


2 


SDSS J143930.77+033317.3 


19.38 


1.80 


1.09 


0.69 


0.58 


0.57 


0.40 


-0.85 


3 


SDSS J145447.32+01 1006.8 


20.35 


1.44 


1.27 


0.74 


0.57 


0.63 


0.55 


-0.57 


3 


SDSS J145547.00+602837.3 


19.18 


1.70 


1.06 


0.63 


0.62 


0.70 


0.47 


-0.28 


2 


SDSS J15051 1.33+620926.3 


18.61 


1.81 


1.07 


0.57 


0.77 


0.57 


0.34 


-0.38 


3 


SDSS J161348.84+482016.0 


18.30 


1.72 


0.96 


0.51 


0.78 


0.64 


0.43 


-0.34 


1 


SDSS J171745.22+625337.0 


18.66 


1.65 


1.24 


0.68 


0.51 


0.68 


0.42 


0.10 


3 


SDSS J173452.52+603603.1 


18.78 


1.71 


1.06 


0.59 


0.66 


0.74 


0.55 


-0.28 


2 


CHCC TT1 CAO-7 tZCi i AACC3£ o 

aUiij) J21593 /. 69+005536. 2 


lo.96 


1 /; 1 
1.61 


1 AO 

1.03 


A C O 

0.5o 


A ^1 

0.63 


A -7A 
0. /O 


A <0 

0.53 




2 


SDSS J221500.88+005217.2 


19.08 


1.60 


1.09 


0.64 


0.51 


0.69 


0.49 


-0.28 


2 


SDSS J221625.03-003122.5 


19.28 


1.58 


1.18 


0.61 


0.57 


0.71 


0.53 


-0.49 


3 


SDSS J223802.82-082532.4 


19.71 


1.78 


1.16 


0.65 


0.49 


0.68 


0.45 


-0.65 


3 


SDSS J224605.41+141640.6 


17.00 


1.57 


1.26 


0.65 


0.49 


0.72 


0.45 


-0.41 


3 


SDSS J224854.83-091723.2 


19.85 


1.57 


1.34 


0.48 


0.54 


0.68 


0.50 


-0.21 


3 


SDSS J225538.34-005945.2 


19.89 


1.47 


0.99 


0.51 


0.61 


0.67 


0.55 


0.21 


1 


SDSS J230303.49-010656.7 


18.95 


1.62 


1.17 


0.63 


0.55 


0.70 


0.50 


-0.40 


3 


SDSS J230805.24+001812.7 


19.36 


1.70 


1.30 


0.72 


0.50 


0.68 


0.42 


-0.49 


3 


SDSS J233030.19+004521.9 


19.45 


1.58 


1.23 


0.58 


0.51 


0.70 


0.46 


-0.80 


3 


SDSS J235 116.25-003104.8 


19.48 


1.49 


1.15 


0.67 


0.53 


0.69 


0.45 


-1.30 


3 


SDSS J235830.60-01 1413.2 


19.99 


1.54 


1.19 


0.60 


0.66 


0.57 


0.49 


-0.61 


2 



a Photometry is omitted when uncertainties are larger than 0.2 magnitudes. 
b Photometry is not available for some objects. 

c EWs of subdwarfs that were either classified as active or inactive (see Section 3.2) are included. 



